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Abstract 

Volatility of intra-day stock market indices computed at various time horizons 
exhibits a scaling behaviour that differs from what would be expected from 
fractional Brownian motion (fBm). We investigate this anomalous scaling 
by using empirical mode decomposition (EMD), a method which separates 
time series into a set of cyclical components at different time-scales. By 
applying the EMD to fBm, we retrieve a scaling law that relates the variance 
of the components to a power law of the oscillating period. In contrast, when 
analysing 22 different stock market indices, we observe deviations from the 
fBm and Brownian motion scaling behaviour. We discuss and quantify these 
deviations, associating them to the characteristics of hnancial markets, with 
larger deviations corresponding to less developed markets. 

Keywords: Empirical mode decomposition, Hurst exponent. Multi-scaling, 
Market efficiency. 


1. Introduction 

Over the last few years hnancial markets have witnessed the availability 
and widespread use of data sampled at high frequencies. The study of these 
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data allows to identify the intra-day structure of financial markets BB- 
Data at these frequencies have dynamic properties which are not generated 
by a single process but by several components that are superimposed on top 
of each other. These components are not immediately apparent, but once 
identihed, they can be meaningfully categorized as noise, cycles at different 
time-scales and trends [l|- 

Since the early work of Mandelbrot B Q > recognized that different 

time-scales contribute to the complexity of hnancial time series in a self¬ 
similar (fractal) manner. Empirical properties of hnancial data at various 
frequencies have been observed in a number of studies, see for example (a, la, 

01 , 0 . 

According to the random walk hypothesis hnancial market dynam¬ 
ics can be described by a random walk, a self-similar process with sca ling 
exponent (Hurst exponent) H = 0.5 [ll|. Opposing this theory, Peters [l2| 
introduced the fractal market hypothesis (FMH), which represents hnancial 
market dynamics by fractional Brownian motion (fBm), a self-similar pro¬ 
cess with scaling exponent Q < H < 1. The focus of the FMH is on the 
interaction of agents with various investment horizons and dihering interpre¬ 
tations of information. Based on this theory, heterogeneous markets models 
have explained some stylised facts (such as volatility clustering, kurtosis, fat 
tails of returns, power law behaviours) observed in hnancial markets, see for 
example [l^, 14, EQ. 


In self-similar uni-scaling process, such as fBm, all time-scales contribute 
proportionally and there is a _specihc relation that links statistical proper¬ 
ties at diherent time-scales [l7|. However, real hnancial time series have 
more complex scaling patterns, with some time-scales contributing dispropor- 
tionally; these patterns characterize multi-scaling processes whose statistical 


properties vary at each time-scale [18|, |19 


21 


3. 


The knowledge of scaling laws in hnancial data helps us to understand 
market dynamics (3 , 0) that can be interpreted to construct efficient and 
prohtable trading strategies. In this paper, we use emp irical mode decompo¬ 
sition (EMD), an algorithm introduced by Huang [25|, to decompose intra¬ 
day hnancial time series into a trend and a hnite set of simple oscillations. 
These oscillations, called intrinsic mode functions (IMFs), are associated with 
the time-scale of cycles latent in the time series. The EMD provides a tool 
for an exploratory analysis that takes into account both the hne and coarse 
structure of the data. This decomposition has been widely used in many 
helds, including the analysis of hnancial time series 2a, l27|, l2a, |29|, , river 
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flow fluctuations [Mi, wind speed [32i, heart rate variability [33l, etc. 

In this paper, we hrst apply EMD to fBm, uncovering a power law scaling 
between the period and variance of the IMFs with scaling exponent related 
to the Hurst exponent. We then apply EMD to 22 different stock market 
indices whose prices are sampled at 30 second intervals over a time span of 6 
months. In this case, we encounter more complex scaling laws than in fBm. 
The deviations from the fBm behaviour are quantihed and interpreted as an 
anomalous multi-scaling behaviour. 

This paper is organized as follows. In section 2, we introduce the EMD. 
In section 3, we present the variance scaling properties of fBm. In section 
4, we present an application to high frequency hnancial data. We hnally 
conclude in section 5. 


2. Empirical mode decomposition 


The empirical mode decomposition is a fully data-driven decomposition 
that can be applied to non-stationary and non-linear data j^. Differently 
from the Fourier and the wavelet transform, the EMD does not require any 
a priori hlter function 3^. The purpose of the method is to identify a hnite 


set of oscillations with scale dehned by the local maxima and minima of the 
data itself. Each oscillation is empirically derived from data and is referred 
as an intrinsic mode function. An IMF must satishes two criteria: 


1. The number of extrema and the number of zero crossings must either 
be equal or differ at most by one. 

2. At any point, the mean value of the envelope dehned by the local 
maxima and the envelope dehned by the local minima is zero. 

The IMFs are obtained through a process that makes use of local extrema 
to separate oscillations starting with the highest frequency. Hence, given a 
time series x{t), t = 1, 2,..., T, the process decomposes it into a hnite number 
of intrinsic mode functions denoted as IMFkit), k = l,...,n and a residue 
Tnit). If the decomposed data consist of uniform scales in the frequency 
space, the EMD acts as a dyadic hlter and the total number of IMFs is 
close to n= log 2 (T) [^. The residue is the non-oscillating drift of the data. 
At the end of the decomposition process, the original time series can be 
reconstructed as: 


x{t) = ^ IMFkit) + Tnit). (1) 

k=l 
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The EMD comprises the following steps: 

1. Initialize the residue to the original time series ro{t) = x{t) and set the 
IMF index k = 1. 

2. Extract the IMF: 

(a) initialize ho(t) = rk-i{t) and the iteration counter i = 1; 

(b) hnd the local maxima and local minima of 

(c) create the upper envelope Eu(t) by interpolating between the max¬ 
ima (lower envelope Ei{t) for minima, respectively); 

(d) calculate the mean of both envelopes as 

(e) subtract the envelope mean from the input time series, obtaining 
hi(t) = hi-i{t) - 

(f) verify if hi{t) satishes the IMF’s conditions: 

• if hi{t) does not satisfy the IMF’s conditions, increase i = i + 1 
and repeat the sifting process from step (b); 

• if hi{t) satishes the IMF’s conditions, set IMEk{t) = hi and 
dehne rfc(t) = rk-i{t) - IMEk{t). 

3. When the residue rk{t) is either a constant, a monotonic slope or con¬ 
tains only one extrema stop the process, otherwise continue the decom¬ 
position from step 2 setting k = k + 1. 


Orthogonality cannot be theoretically guaranteed, but in most cases it 
is satished j^. Including the residue as the last component and rewriting 

n+1 

Equation [T] as x{t) = Ckit), the square of the values of x{t) can be 


expressed as: 


k=l 


n+1 


n+l n+l 


k=l j=l k=l 

If the decomposition is orthogonal, the cross terms should be zero. An 


index of orthogonality (10) is dehned as [25 


io=j: 


n+l n+l 

E E c,(,t)c,{t) 

j=l k=l 


x{ty 


(3) 
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3. Self-similar scaling exponent 


Self-similarity or scale invariance is an attribute of many laws of nature 
and is the underlying concept of fractals. Self-similarity is related to the 
occurrence of similar patterns at different time-scales. In this sense, proba¬ 
bilistic properties of self-similar processes remain invariant when the process 
is viewed at different time-scales j^, 37, 3^ . 

A stochastic process X{t) is statistically self-similar, with scaling expo¬ 
nent 0 < 77 < 1, if for any real a > 0 it follows the scaling law: 


X{at) = a^X{t) 


t e 


(4) 


where the equality (=) is in probability distribution 38 


An example of self-similar process is fractional Brownian motion (fBm), 
a stochastic process characterized by a positive scaling exponent 0 < 77 < 1 
39|. When 0 < 77 < h fBm is said to be anti-persistent with negatively 


auto-correlated increments. For the case 4 < 77 < 1, fBm reflects a persistent 
behaviour and the increments are positively auto-correlated. When 77 = 4, 
fBm is reduce to a process with independent increments known as Brownian 
motion. 


3.1. EMD based scaling exponent 

Flandrin et al. ji^ empirically showed that when decomposing fractional 
Gaussian noise (fGn), the differentiation process of fBm j^, the EMD can be 
used to estimate the scaling exponent 77, if 77 >i The authors ascertained 
that the variance progression across IMFs satishes, var{IMF^^"') oc 
where the function denotes the period of the 7*^-IMF 0. 

In this paper, we follow a similar approach to ji^ , but instead of applying 
EMD to fGn, we considered its integrated process, fBm. We empirically 
showed that a similar scaling law holds for the variance of IMFs: 

var{IMFl^'^) oc . (5) 

Therefore, for fBm and fGn, the IMF variance follows a power law scaling 
behaviour with respect to its particular period of oscillation, and the scaling 


'^The periods can be approximated as the total number of data points divided by 
the total number of zero crossings of each IMF. 
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parameter is related to the Hurst exponent. The EMD estimator of H can 
be determined by the slope of a linear regression £t on the logarithmic of the 
variance as a function of the logarithmic of the period, 


log log(Ti) + log(co). 


( 6 ) 


where Cq is the intercept constant of the linear regression. In the following 
section we will provide the simulation results supporting equation |6l 

3.2. FBm simulation analysis 

In order to verify the empirical scaling law of Equation El we generated 
N = 100 fBm processes for the following values of the scaling exponent 
H = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9. The simulated processes have 
two different lengths, Ti = 10, 000 and T 2 = 100, OOoJl 

We applied the EMD to each fBm simulation and calculated its respec¬ 
tive H* exponent. In Table [221 we report the mean over the 100 

estimators. We also report the root mean square error (RMSE) of the estima¬ 



tors, RMSE = y ——-. We observe that the longer the analysed time 

series, the better the estimation of H is. For length T = 100, 000, is 

indeed very close to the scaling exponent H (for all values of H). In Figure 
in we plot the mean values of the H* exponent as presented in Table 13.21 
The error bars represent the RMSE of the estimator. Let us emphasize that 
we do not propose the EMD as a way to estimate the Hurst exponent, but 
as a tool to analyse the interactions between the different time scales present 
in the data. For comparison, we estimated the Hurst exponent using the 
generalized exponent approach jl^ . In Table 13.21 we include the mean and 
the RMSE of this estimator denoted as He- 

Moreover, to visualize the linear relationship of Equation [6l we explicitly 
show the relation between log(nar(/MF/'®™')) and log(rfc) for a fBm simula¬ 
tion of scaling exponent H = 0.6 and length Ti = 10, 000 points, see Figure 
|2] . In this example, the resulting estimator is H* = 0.593 which accurately 
approximates the scaling exponent of the simulated process. 


^All the fflm paths were generated using MATLAB® wavelet toolbox. 
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10,000 100,000 


H 

{H*) 

RMSEh^ 

{Hg) 

RMSEho 

( ff *) 

RMSEh- 

(Hg) 

RMSEua 

0.1 

0.05 

0.06 

0.15 

0.05 

0.11 

0.03 

0.15 

0.05 

0.2 

0.15 

0.07 

0.22 

0.03 

0.21 

0.03 

0.22 

0.03 

0.3 

0.26 

0.06 

0.31 

0.01 

0.31 

0.03 

0.31 

0.01 

0.4 

0.38 

0.05 

0.40 

0.01 

0.41 

0.02 

0.40 

0.01 

0.5 

0.49 

0.05 

0.50 

0.01 

0.51 

0.03 

0.50 

0.01 

0.6 

0.59 

0.04 

0.60 

0.01 

0.60 

0.03 

0.60 

0.01 

0.7 

0.70 

0.04 

0.70 

0.01 

0.69 

0.03 

0.70 

0.01 

0.8 

0.80 

0.04 

0.79 

0.01 

0.78 

0.04 

0.79 

0.02 

0.9 

0.90 

0.05 

0.88 

0.03 

0.87 

0.04 

0.87 

0.03 


Table 1: Confirmation that the empirical scaling law of Eq. 3 retrieves the expected 
scaling exponent for fractional Brownian motion. Mean of the scaling exponent H* over 
100 simulations of fBm with parameter H = 0.1,0.2,..., 0.9 and length; left Ti = 10,000 
and right: T 2 = 100, 000. For comparison, we included the mean and the RMSE of the 
generalized Hurst exponent estimator denoted as Hq- 




Figure 1: Demonstration that the empirical scaling law of Eq. 3 retrieves the expected scaling 
exponent for fractional Brownian motion. Mean of the scaling exponent H* over 100 simulations 
of fBm with parameter H = 0.1,0.2,..., 0.9 and length: left Ti = 10,000 and right: T 2 = 100,000. 
The error bars denote the RMSE of the estimator. 
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Figure 2: Log-log plot of IMF variance as a function of period for a fBm oi H = 0.6 and 
length Ti = 10,000. The blue line represents the least square fit. The scaling exponent 
H* = 0.593 can be recovered from half the slope of the least square fit. 


4. Variance scaling in intra-day financial data 

We analysed intra-day prices for 22 different stocks market indices. The 
data set, extracted from Bloomberg, covers a period of 6 months from Mw 
5 */!, 2014 to November 5*^, 2014. Prices are recorded at 30 second interval^. 
We exclnded weekends and holidays. The nnmber of working days and the 
nnmber of points for every trading day depend on the opening honrs of each 
stock exchange. The list of analysed stock market indices is reported in Table 

El 

We applied EMD to the logarithmic price of each financial time series. 
For the sake of clarity, in this section we only focus on the decomposition of 
the S&P 500 index, but a similar analysis has been done for the other stock 
market indices. For the S&P 500 log-price time series, we extracted 17 IMFs 
and a residue that describe the local cyclical variability of the original signal 
and represent it at different time-scales. The original log-price time series 


^Prices for the Warsaw stock exchange index (WIG) were only available at every minute 
frequency. 




Country 

Index 

Length 

Brazil 

BOVESPA 

105,000 

China 

SSE 

60,480 

France 

CAC 40 

136,080 

Greece 

ASE 

106,470 

Hong Kong 

HSI 

98,154 

Hungary 

BUX 

122,880 

Italy 

FUSE MIB 

133,056 

Japan 

NIKKEI 225 

75,600 

Malaysia 

KLSE 

115,320 

Mexico 

IPC 

100,620 

Netherlands 

AEX 

130,680 

Poland 

WIG 

64,680 

Qatar 

DSM 

52,080 

Russia 

RTSI 

133,120 

Singapore 

STI 

123,840 

South Africa 

JSE 

117,500 

Spain 

IBEX 

135,527 

Turkey 

XU 100 

91,760 

UAE 

UAED 

60,000 

UK 

FUSE 

130,560 

USA 

S&P 500 

99,840 

USA 

NASDAQ 

100,620 


Table 2: Stock market indices including the length of the time series. 


and its IMFs are displayed in FignreEi In this figure, we observe temporary 
clusters of volatility that characterize some of the components, for example 
the high volatility at the end of the time series can evidently be seen in 
components 2,3,4,6 and 7. 

The IMF periods, calculated as the total number of data points divided by 
the total number of zero crossings, are reported in Table O These periods are 
converted into minutes, hours and days. The fastest component has a cycle 
of 1.6 minutes, contrasting the slowest cycle of 11.6 days. Notice that the 
first 12 IMFs represent the intra-day activity (6.5 hours of trading), while the 
remaining IMFs (from 13*^ to 17*^) are associated with the inter-day cycles. 
The last component is the residue of the EMD. 

The overall trend of the time series is given by the residue, and each 
component can be seen as an oscillating trend of the previous component 
on a shorter time-scale. The effectiveness of EMD as a de-trending and 
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Figure 3: Top: log-price time series of the S&P 500 index for the period 05/05/2014 
to 05/11/2014. Bottom: the 17 IMFs and the residue obtained through EMD of the 
log-prices. 
















































































































IMF 

Period/min 

IMF 

Period/hr 

IMF 

Period/days 

1 

1.6 

9 

1.1 

14 

1.1 

2 

2.8 

10 

1.9 

15 

2.2 

3 

4.9 

11 

3.0 

16 

4.3 

4 

8.4 

12 

5.9 

17 

11.6 

5 

13.0 

13 

11.7 

18 

Residue 

6 

19.3 





7 

28.8 





8 

41.7 






Table 3: Period of the IMFs obtained from the S&P 500 index. 



Figure 4: Log-price time series of the S&P 500 index (blue line). The red line represents 
the ‘trend’ of the data calculated as the sum of the residue plus the last IMF. 


smoothing tool is illustrated in Figure 01 In this hgure, the original time 
series (blue line) is compared with a ’trend’ (red line), calculated as the sum 
of the residue plus the last component. 

In previous section we discussed that for fBm, the EMD produces a lin¬ 
ear relationship between the logarithmic values of variance of the IMFs and 
its respective period of oscillation (Equation 0]). We tested whether this re¬ 
lationship also holds for hnancial data. In Figure 01 we show the log-log 
plot of variance as a function of period for the IMFs obtained from the S&P 
500 index (red diamonds). The estimated scaling exponent has a value of 
H* = 0.55. The goodness of the linear £t was estimated by the coefficient of 
determination!^ which is = 0.992. We can conclude that this index satishes 


^This coefficient of determination is the square of correlation between the dependent 
and independent variable. Values of this coefficient range from 0 to L with 1 indicating a 
perfect ht between the data and the linear model, see for example |41j . 
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Figure 5: Log-log plot of IMF variance as a function of period for the EMD of the S&P 
500 index. The red line represents the best least square linear fit. The goodness of the 
linear fit is = 0.992. 


the linear relationship of Eqnation|6]bnt the scaling exponent, H* = 0.55, is 
different from that of Brownian motion, H = 0.5. 

We performed the same analysis for the other stock indices. Ending both 
significant deviations from Brownian motion [H* ^ 0.5 ) and deviations from 
the scaling law of Eqnation El 

In Table 01 we report more details abont the decomposition for each fi¬ 
nancial index. We inclnde the nnmber of IMFs and the index of orthogonality 
as described in Eqnation El We observe small valnes of the 10, indicating 
an almost orthogonal decomposition. Fnrthermore, we report the estimated 
exponent and the goodness of the linear fit for every stock market index. 
Althongh the coefficients of determination are all above 0.94, we shall dis- 
cnss shortly that significant deviations from linearity (fBm behavionr) are 
observed, especially in less developed markets. 

Let ns now discnss in more detail the deviations of the scaling laws fonnd 
in stock markets from the scaling expected in Brownian motion (Bm). With 
this aim, we generated N = 100 paths of Bm with length T eqnal to the 
analysed stock market index (see Table [2]). We applied EMD to each simn- 
lation and obtained its respective intrinsic oscillations denoted as IMF^^\ 
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Index 

# IMFs 

lO xl0'‘ 

W 

H* 

S&P 500 

17 

3.7 

0.992 

0.564 

BOVESPA 

18 

6.3 

0.989 

0.561 

FTSE MIB 

18 

8.6 

0.987 

0.571 

XU 100 

19 

3.5 

0.985 

0.563 

RTS 

20 

11 

0.985 

0.581 

CAC 40 

17 

00 

00 

0.984 

0.564 

UAED 

16 

6.3 

0.978 

0.616 

ETSE 

23 

2.9 

0.977 

0.529 

ASE 

15 

29 

0.974 

0.587 

IBEX 

18 

5.7 

0.973 

0.531 

WIG 

16 

8.2 

0.973 

0.591 

SSE 

14 

1.6 

0.971 

0.534 

DSM 

18 

7.6 

0.971 

0.618 

IPG 

18 

0.68 

0.971 

0.555 

BUX 

19 

4.2 

0.970 

0.542 

HSI 

19 

1.0 

0.969 

0.554 

AEX 

21 

13 

0.968 

0.558 

NASDAQ 

20 

5.1 

0.960 

0.530 

NIKKEI 225 

22 

8.4 

0.959 

0.544 

JSE 

19 

2.9 

0.956 

0.518 

KLSE 

22 

0.77 

0.943 

0.540 

STI 

21 

2.6 

0.942 

0.522 


Table 4: Stock market indices including the number of IMFs obtained when applying EMD 
to the logarithmic price. The second column report the index of orthogonality (xlO"^). 
Stock market indices are reported in descending order of which represents the goodness 
of the linear fit of Equation [H] Last column reports the estimated exponent H* of the 
same equation. 


i = 1, 2,..., 100, k = 1, 2,..., rtj, with rii the number of IMFs for each Bm 
simulation. In order to compare the variance of the IMFs extracted from the 
hnancial index X, var{IMF^), against the var^IMF^'^^)^ we rescaled the 
latter as: 


var{IMF^^^) = a var{IMF^^^), 

where 


n 

iT.{var(IMFi^)/T^) 

k=l 

rii 

k=l 


(7) 

( 8 ) 
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Figure 6: Log-log plot of variance as a function of period for the S&P 500 IMFs (red 
diamonds) compared with 100 rescaled Bm linear fits of slope H* = 0.5 (blue lines). 


For the rescaled var{IMFj^^'), we estimated the intercept constant Co^ 
of Eqnation [6l hxing H = 0.5. In Fignre El we present all these 100 linear 
hts as light blue lines. In the same hgure, we plotted the variance of the 
IMFs extracted from the S&P 500 index, same as reported in Figure El We 
observe that the Brownian motion linear hts (blue lines) and the linear £t of 
the S&P 500 index (red line) are close to each other, suggesting an efficient 
behaviour in this market. 

The goodness of the linear £t between the hnancial data points (red dia¬ 
monds) and each of the Brownian motion linear £t (blue lines) was calculated 
as follow: 


n r, 

^ [log {var (IMF^)) - log (q co,rf)] 

Rim. = 1 - --• ( 9 ) 

E [log (var (IMF^)) - (log (var (IMF^^)))]" 


The number of IMFs extracted from the stock index X is denoted by n and 
(■) indicates the mean over these n IMF variances. The deviations from 
Bm were calculated by the mean over the goodness of the linear hts, i.e. 
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100 

we calculated {R\rn) = “ X] ^Bm ■ For the S&P 500 index, we obtained a 

i=\ 

coefficient equal to (-R^m) = 0.979, demonstrating the similarity between the 
scaling properties of this hnancial index and Brownian motion. 


4.1- Scaling properties of stock market indices 

In previous section, we introduced two measures that quantify the devia¬ 
tions from the scaling behaviour of fractional Brownian motion and Brownian 
motion. These measures are given by: 

1. coefficient of determination (square of correlation) between the 
logarithm of period and logarithm of variance of IMFs obtained from 
stock market indices. 

2. mean of the relative squared residuals between the IMF vari¬ 
ances obtained from hnancial data and each of the linear hts for Brow¬ 
nian motion simulations. 


In Table [5l we report the values of ^Fe 22 stock market indices. 

For comparison purposes, we repeated the values. The last column in this 
table indicates the ordering of the markets if R? were used as the ranking 
measure. 

The S&P 500 index is ranked the highest in both scales. Developed mar¬ 
kets tend to be at the top of the table with some exceptions that may arise 
from the specihc characteristics of the analysed period of time. May 5*^, 
2014 to November 5*^, 2014. In Figure [TJ we plot the hnancial market rank- 


The horizontal bars represent the 5*^ and 95*^ percentiles of the R 


ing. J-liC liUi IZ/UlibCXl uaio iCpiCDCllL bliC O aiiU. ao pCi L-ClibllCD Ui bliC 

distribution. The blue dot inside each bar indicates the mean value {R\rn) 
as reported in Table Ei Despite the fact that some hnancial stock indices 
have similar values of R^rm^ recognize statistically signihcant diher- 

ences between developed and emerging markets, observing a clear tendency 
for the developed markets to present larger values of {R%^ with narrower 
distributions. 

In order to visualize the anomalous scaling in some stock markets and to 
understand the origin of the diherences in the results, we compare the cases 
of the NASDAQ (USA), BOVESPA (Brazil), NIKKEI 225 (Japan) and DSM 
(Qatar) indices in more detail. 

For the NASDAQ index (USA), we obtained 20 IMFs and a residue. In 
Figure 8(a), we present the log-price time series (blue line) and the ‘trend’ 
consisting of the residue plus the last IMF (red line). In Figure 8(b), we 
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Country 

Index 


Rank^j^2^ ^ 


Rankp2 

USA 

S&P 500 

0.979 

1 

0.992 

1 

Brazil 

BOVESPA 

0.977 

2 

0.989 

2 

UK 

FTSE 

0.973 

3 

0.977 

8 

Turkey 

XU 100 

0.972 

4 

0.985 

4 

Italy 

FTSE MIB 

0.971 

5 

0.987 

3 

France 

CAC 40 

0.970 

6 

0.984 

6 

Spain 

IBEX 

0.969 

7 

0.973 

10 

China 

SSE 

0.967 

8 

0.971 

12 

Russia 

RTSI 

0.964 

9 

0.985 

5 

Hungary 

BUX 

0.963 

10 

0.970 

15 

Mexico 

IPC 

0.960 

11 

0.971 

14 

Hong Kong 

HSI 

0.958 

12 

0.969 

16 

USA 

NASDAQ 

0.954 

13 

0.960 

18 

Netherlands 

AEX 

0.953 

14 

0.968 

17 

South Africa 

JSE 

0.952 

15 

0.956 

20 

Japan 

NIKKEI 225 

0.949 

16 

0.959 

19 

Greece 

ASE 

0.948 

17 

0.974 

9 

Poland 

WIG 

0.947 

18 

0.973 

11 

UAE 

UAED 

0.939 

19 

0.978 

7 

Singapore 

STI 

0.934 

20 

0.942 

22 

Malaysia 

KLSE 

0.933 

21 

0.943 

21 

Qatar 

DSM 

0.928 

22 

0.971 

13 


Table 5: Stock market indices ranked in descending order of {R%rn)- The last column 
indicates the ordering of the markets with respect to . 
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DSM 
KLSE 
STI 
UAED 
WIG 
ASE 
NIKKEI 225 
JSE 
AEX 
NASDAQ 
HSI 
I PC 
BUX 
RTS I 
SSE 
IBEX 
CAC 40 
FTSE MIB 
XU 100 
FTSE 
BOVESPA 
S&P 500 

0.92 0.93 0.94 0.95 0.96 0.97 0.98 

Figure 7: Percentiles 5*^ and 95*^ of the distribution for the analysed stock market 

indices. The blue dot inside each bar indicates the value of {R^sm) used for the financial 
market ranking. 
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(a) (b) 


Figure 8: EMD analysis for the NASDAQ index. Captions for figures (a) and (b) are the 
same as captions for figures (4) and (6) respectively. 


observe that the deviation from the linear relationship of Equation [6] is sig- 
nihcant. Thus, the log-log relationship between period and variance is not 
completely satished. The resultant coefficient of determination is = 0.960, 
ranking this index at the 18*^ position. Moreover, when compared with Bm, 
we identify that most of the components deviate from the Bm linear hts (blue 
lines). We also note that the total number of component (21) is considerable 
larger than what would be expected from a process with uniform scales, i.e., 
Iog2(100620)=16.6. The presence of these extra oscillations with reduced 
variance suggests a more complex structure than Bm. The deviations from 
Bm model, quantihed by the coefficient {Rsm) — 0.958, rank this index at 
the 13^^ position. 

The variance scaling properties of the BOVESPA index (Brazil) are pre¬ 
sented in Figure ini For this stock index, the EMD identihes long period cycles 


with larger variance than what would be expected from Bm, see Figure 9(b) 


However, the linear £t between the logarithmic value of IMF variances and 
periods is in general good with R^ = 0.989. The goodness of the linear fit 
between the Bm simulations is {Rsm) — 0.977, placing this index at the 
second position. Such a good ranking for this market may be unexpected, 
but we must stress that it only reflects the six-month period of observations. 
From Figure 9(a), we can see that this was a rather random but calm period. 
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(a) (b) 


Figure 9: EMD analysis for the BOVESPA index. Captions for figures (a) and (b) are 
the same as captions for figures (4) and (6) respectively. 


For the NIKKEI 225 index (Japan), we obtained 22 IMFs and a residue. 
Similar as the NASDAQ index, the number of components is considerable 
larger than what would be expected from Bm, i.e., Iog2(75600)=16.2. These 
many oscillations, specially the high frequency components, generate a non¬ 
linear behaviour that deviates from Bm. Given the anomalous scaling be¬ 
haviour of this stock index, see Figure 10(b), we obtained {Rsm) — 0.949, 
ranking it at the 16*^ position. 

Finally, the DSM stock index (Qatar) is displayed in Figure dH The log- 
price time series and its respective ‘trend’are displayed in Figure 11(a) In 


Figure 11(b), we observe the poor liner £t of Equation [6] that is characterized 
by a considerable steep slope. We obtained = 0.971, ranking this index 
at the lowest position. Furthermore, if we compare its IMF variances against 
the Bm linear hts, we observe that most of the variance values (red diamonds) 
follow outside the band expected from Bm. The large variance of the low 
frequency components suggests the presence of important long period cycles. 
Given its deviations from Bm behaviour, this index is also ranked the lowest 
with respect to the measure = 0.928. 


5. Conclusions 

We explored the scaling properties of EMD, an algorithm that detrends 
and separates time series into a set of oscillating components called IMFs 
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(a) (b) 

Figure 10: EMD analysis for the NIKKEI 225 index. Captions for figures (a) and (b) 
are the same as captions for figures (4) and (6) respectively. 



(a) (b) 


Figure 11: EMD analysis for the DSM index. Captions for figures (a) and (b) are the 
same as captions for figures (4) and (6) respectively. 
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which are associated with specific time-scales. We empirically showed that 
fBm obeys a scaling law that relates linearly the logarithm of the variance 
and the logarithm of the period of the IMFs. For fBm, we demonstrated 
that the extracted coefficient of proportionality equals the scaling exponent 
H multiplied by two. When applied to stock market indices, the EMD re¬ 
veals instead different scaling laws that can deviate significantly from both 
Brownian motion and fractional Brownian motion behaviour. In particular, 
we noted that the EMD of high frequency financial data results in a larger 
number of IMFs than what would be expected from Brownian motion. These 
many components, specially with high frequencies, create a curvature that 
disobeys the linearity in the log-log relation between IMF variance and period 
found in fBm. This is a direct indication of anomalous scaling that reveals a 
more complex structure in financial data than in self-similar processes. 

In this study, we applied EMD to 22 different stock indices and observed 
that developed markets (European and North American markets) tend to 
have scaling properties closer to Brownian motion properties. Conversely, 
larger deviations from uni-scaling laws are observed in some emerging mar¬ 
kets such as Malaysian and Qatari. 

These findings are in agreement with the discernible characteristics of 
developed and emerging markets, the former t ype being more likely to ex¬ 
hibit an efficient behaviour, see for example [l8|, l42| . Compared to previous 
approaches, the EMD method has the advantage to directly quantify the 
cyclical components with strong deviations, giving a further instrument to 
understand the origin of market inefficiencies. 
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